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21 ABSTRACT 

22 The objective of this study was to compare several approaches to soil moisture (SM) retrieval 

23 using L-band microwave radiometry. The comparison was based on a brightness temperature 

24 (Tb) data set acquired since 2010 by the L-band radiometer ELBARA-II over a vineyard field 

25 at the Valencia Anchor Station (VAS) site. ELBARA-II, provided by the European Space 

26 Agency (ESA) within the scientific program of the SMOS (Soil Moisture and Ocean Salinity) 

27 mission, measures multiangular Tb data at horizontal and vertical polarization for a range of 

28 incidence angles (30°-60°). Based on a three year data set (2010-2012), several SM retrieval 

29 approaches developed for spacebome missions including AMSR-E (Advanced Microwave 

30 Scanning Radiometer for EOS), SMAP (Soil Moisture Active Passive) and SMOS were 

31 compared. The approaches include: the Single Channel Algorithm (SCA) for horizontal 

32 (SCA-H) and vertical (SCA-V) polarizations, the Dual Channel Algorithm (DCA), the Land 
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Parameter Retrieval Model (LPRM) and two simplified approaches based on statistical 
regressions (referred to as 'Mattar' and 'Saleh'). Time series of vegetation indices required for 
three of the algorithms (SCA-H, SCA-V and ‘Mattar’) were obtained from MODIS 
observations. The SM retrievals were evaluated against reference SM values estimated from a 
multiangular 2-Parameter inversion approach. The results obtained with the current base line 
algorithms developed for SMAP (SCA-H and -V) are in very good agreement with the 
‘reference’ SM data set derived from the multi-angular observations (R 1 2 ~ 0.90, 
RMSE varying between 0.035 and 0.056 m 3 /m 3 for several retrieval configurations). This 
result showed that, provided the relationship between vegetation optical depth and a remotely- 
sensed vegetation index can be calibrated, the SCA algorithms can provide results very close 
to those obtained from multi-angular observations in this study area. The approaches based on 
statistical regressions provided similar results and the best accuracy was obtained with the 
‘Saleh’ methods based on either bi-angular or bipolarization observations (R 2 = 0.93, 
RMSE ~ 0.035 m 3 /m 3 ). The LPRM and DCA algorithms were found to be slightly less 
successful in retrieving the 'reference' SM time series (R 2 ~0.75, RMSE ~ 0.055 m 3 /m 3 ). 
However, the two above approaches have the great advantage of not requiring any model 
calibrations previous to the SM retrievals. 


1. Introduction 

Surface soil moisture plays a major role in the water and energy budgets of continental 
surfaces, which has direct implications for hydrological, climate, and weather forecasting 
models. L-band passive microwave remote sensing is one of the most promising approaches 
to monitor this variable at the global scale with frequent revisiting times (Jackson et al., 1995; 
Kerr et al., 2001, Njoku et al., 2003; De Lannoy et al., 2013). Three recent or planned space 
missions use this technology: SMOS (launched end of 2009), Aquarius (launched in June of 
2011) and SMAP (launch scheduled in November 2014). 

The Soil Moisture and Ocean Salinity (SMOS) mission is the first spaceborne mission 

dedicated to soil moisture (SM) mapping. SMOS has multi-angular capabilities which are 

exploited by the SM retrieval approach: SM and vegetation optical depth x (used to 

parameterize vegetation attenuation and emission) are retrieved simultaneously based on 

SMOS multi-configuration observations, in terms of polarizations and incidence angles. 

Aquarius is a combined passive/active L-band microwave instrument which consists of a set 
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of three radiometers and a scatterometer, operating at 1 .4 GHz and 1 .26 GHz respectively 
(Levine et al., 2010). The primary mission objective of Aquarius is to provide global 
observations of surface sea salinity once every 7 days. However, Aquarius has also potential 
capabilities to monitor soil moisture at global scales (Luo et al., 2013, Bindlish et al., 2013). 
SMAP incorporates a radar and a radiometer, both operating at L-band and at the incidence 
(observation) angle 0 = 40°. The spatial resolutions of the corresponding active- and passive 
microwave signatures are ~ 39 km x 47 km and ~ 1 km x 1 km, respectively. The mission 
concept is to combine the complementary attributes of the radar observations (high spatial 
resolution but lower soil moisture accuracy) and radiometer observations (higher soil moisture 
accuracy but coarse spatial resolution) to retrieve SM at a spatial resolution of 9 km, and the 
freeze-thaw state at a spatial resolution of 3 km (Entekhabi et al., 2010; O’Neill et al., 2013). 
Several SM retrieval approaches have been developed in the context of these L-band space 
missions. As noted above, in the operational SMOS SM retrieval algorithm, SM and 
vegetation optical depth at nadir (tnad) are retrieved simultaneously based on SMOS 
multiangular and bipolarization observations (Wigneron et al., 1995, 2000; Kerr et al., 2012). 
The 2-Parameter (2-P) retrievals of SM and xnad are obtained from inversion of the L-MEB 
(L-band Microwave Emission of the Biosphere) model (Wigneron et al., 2007). This forward 
model is based on the so-called x-co model (Mo et al., 1982) and it includes a number of 
parameterizations to capture effects of vegetation structure and soil roughness on polarization 
and angular properties of L-band Tb emitted from land surfaces. The inversion of L-MEB 
considering SM and xnad as the requested parameters (referred to as 'L-MEB 2-P' inversion) 
is implemented in the operational algorithms used to compute the Level 2 (distributed by 
ESA) and Level 3 (distributed by the Centre Aval de Traitement des Donnees SMOS 
(CATDS), Berthon et al., 2012) SMOS products. In parallel to this operational retrieval 
method, several simplified methods have been developed to exploit the capability of L-band 
radiometers to provide information on land surface states such as SM. Lor instance, Wigneron 
et al. (2004) and Saleh et al. (2006) have evaluated statistical regressions based on bi- 
polarization or bi-angular T B data. Mattar et al. (2012) have evaluated similar regression 
methods that also use a vegetation index estimated from ancillary remotely sensed 
observations (such as the Normalized Difference Vegetation Index (NDVI) or the Leaf Area 
Index (LAI)) to account for vegetation effects. Moreover, methods based on Neural Networks 
have been and are currently evaluated (Liu et al., 2002; Rodriguez et al., 2003). 
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The general retrieval approach proposed for SMAP is different from the operational SMOS 
SM retrieval: SMAP observations will be available for the sole incidence angle of 40°, but 
make use of the complementary information provided by the active- (radar) and the passive 
(radiometer) L-band data. In the initial release of the ATBD (Algorithm Theoretical Basis 
Document) written for the retrievals from SMAP’s radiometer (O’Neill et al., 2013), four soil 
moisture retrieval algorithms are suggested for evaluation during the pre- and post-launch 
calibration and validation activities: (i) the single-channel algorithm at H polarization (SCA- 
H) which is the current SMAP baseline algorithm, (ii) the single-channel algorithm at V 
polarization (SCA-V), (iii) the dual-channel algorithm (DCA), and (iv) the Land parameter 
retrieval model (LPRM). In the SCA-H and -V algorithms, vegetation is accounted for by the 
t-co model as in L-MEB. However, optical depth at nadir (tnad) is not retrieved as for SMOS. 
Instead it is estimated from the linear relation tnad = b ■ VWC between tnad and vegetation 
water content (VWC) (Jackson et al. (1991)). Thereby, values of the b-parameter are assumed 
polarization independent and will be provided from a land cover look up table, and the VWC 
is estimated from values of the NDVI Index. The DCA retrieval approach is very similar to 
the one used for SMOS. The only difference is that the inversion is based on the minimi z ation 
of a cost function accounting for the Root Mean Square Error (RMSE) between measured and 
simulated bi-polarized Tb observations at one incidence angle, whereas multi-angular 
observations are used for SMOS. In the LPRM algorithm, the Microwave Polarization 
Difference Index (MPDI) and the observed emissivities are used to derive the vegetation 
optical depth x (Meesters et al., 2005). In a second step, SM is retrieved with an optimization 
routine that minimizes the error between the modelled and observed H-polarized brightness 
temperatures (Owe et al., 2008; De Jeu et al., 2009). 

In this study, these different retrieval algorithms were compared using a 3 -year long 
multiangular T B data set acquired by the L-band radiometer ELBARA-II over a vineyard field 
(MELBEX-III) at the Valencia Anchor Station (VAS) site (Schwank et al., 2012, Wigneron et 
al., 2012). Applications of the retrieval methods can be made at large scales from satellite 
observations but also at more local scale for long term SM monitoring from ground based 
instruments mounted on different types of platforms: towers as for ELBARA-II (de Rosnay et 
al., 2006; Schwank et al., 2012; Schlenz et al., 2012, etc.); trucks (Hornbuckle et al., 2004; 
Kurum et al., 2009) or from the top of a mountain as in Pellarin et al. (2013). 

ELBARA-II (Schwank et al. 2010), developed by GAMMA Remote Sensing AG 
(Switzerland) and funded by the ESA, provides T B at horizontal and vertical polarization for a 
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range of observation angles (30°-60°). The ELBARA-II Tb observations were acquired since 
2010 and a 3-year Tb data set is available for the MELBEX-III site. As an accurate estimation 
of SM from ground based measurements over the MELBEX-III site could not be achieved 
because of very frequent agricultural practices within the field, it was considered that 
representative SM values (referred to as 'reference' SM data set) over the ELBARA-II 
footprints were obtained from multi-angular 2-P L-MEB retrievals. Moreover, the 2-P L-MEB 
approach also provided retrievals of optical depth at nadir (tnad). These latter values were 
used to calibrate the relationships between tnad and NDVI, which are required in the SCA-H 
and SCA-V algorithms. Based on these ‘reference’ SM and tnad data sets and the ELBARA- 
II T b observations, seven SM retrieval approaches were evaluated and compared: the four 
methods considered presently in the SMAP ATBD based on bi-polarization observations at 
one observation angle (0 = 40° for SMAP) and three regression methods (Saleh et al, 2006 
and Mattar et al., 2012) developed in the framework of SMOS research activities and based 
on bi-angular or bipolarization observations. The results of this evaluation are discussed in the 
context of the improvement and development of the SM retrieval algorithms. 


2. Materials and method 

2.1. The ELABARA-II radiometer at MELBEX-III (VAS site) 

The study was based on Tb measurements made by the ELBARA-II radiometer over the 
2010-2012 period within the VAS site. ELBARA-II was installed in September 2009 at the 
MELBEX-III vineyard field (referred to as M-III), close to Caudete de las Fuentes, on the 
Utiel-Requena Plateau at ~ 800 m a.s.l., in the region of Valencia, Spain (39°31'18.18"N, 
1°17'29.64"W). This site is one of the reference sites selected by ESA in Europe within the 
SMOS science program. 

All details concerning the ELBARA-II instrument and the M-III experiment set up are given 
in Schwank et al. (2010, 2012), and Wigneron et al. (2012). Only a brief summary of the main 
information concerning this experiment is presented here. 

The ELBARA-II radiometer was set up 17 meters above ground to monitor a vineyard that is 
representative of the main land use of the VAS region. The ELBARA-II was equipped with 
an elevation tracker that allows measurements at specific observation angles 0 varying 
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between 30° < 0 < 330° with 0 = 180° being the zenith direction. Every 30 minutes, 
automated "elevation scans" are carried out that provide T B at horizontal and vertical 
polarizations at observation angles between 0 = 30° and 70° with steps of 5°. Between each 
elevation scan, measurements are made at the 0 = 45° every 10 minutes. Once a day, at 23:55 
local time, the radiometer is automatically positioned at 150° to carry out sky calibration 
measurements. The absolute accuracy of the ELBARA-II measurements was estimated to be 
better than ±1 K over the course of 2010-2012. During short time periods, no measurement 
could be acquired over the vineyard field due to experiments using reflecting foils (Schwank 
et al., 2012) or due to technical issues: in 2010 (DoY 222 - DoY 245, DoY 312 - DoY 337) 
and in 2011 (DoY 41 - DoY 62; DoY 84 -DoY 133). The ELBARA-II observations were 
slightly affected by Radio Frequency Interferences (RFI) caused by active microwave systems 
violating the protected part of the L-Band (1400 MHz - 1427 MHz). Efforts made by the 
Spanish administrative authorities in 2010 to mitigate RFI disturbances resulted in a 
significant decrease since the beginning of July in 2010 (~ DoY 190). Most RFI events result 
in steep increases in the time variations of the measured T B (larger than 30K at minimum) and 
unrealistic T B values (larger than 330 K). These RFI events were detected manually from the 
ELBARA-II T b data set. To be consistent with the overpass times of SMOS and SMAP, only 
T b measurements made at 6 am and 6 pm local time are considered in this study. 


2.2. In situ measurements 

Concurrent with the ELBARA-II observations, ground measurements were obtained within 
the M-III vineyard. Soil profiles of the volumetric soil moisture [m 3 m°] and temperature 
were acquired up to about 1 m (Wigneron et al., 2012). Vineyard cultivation practices are 
carried out frequently within the field (for weeding and pest control, winter and summer 
pruning, cluster thinning, etc.) so that SM probes could not be installed permanently within 
the ELBARA-II footprints. Only two Delta-T Theta Probes measuring the volumetric SM of 
the top 0-6cm soil layer were installed at the border of the field where no field work was 
carried out. It is our opinion that these SM probes cannot provide SM values representative of 
the field conditions as seen by the ELBARA-II instrument and have not been used in the 
analysis presented here (Wigneron et al., 2012). 
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A meteorological station located at the VAS (coordinates: 39°34’15”N, 1°17’18”W, 813 m 
a.s.l.), a few kilometres from the M-III site provided the standard meteorological variables (air 
temperature, wind speed, air humidity, etc.)- Over the VAS site, the average value of the total 
yearly precipitation over the ten years prior to 2010 is P = 461 mm. For the three years 
considered in this study; 2010 was wet (P= 538.2 mm) and was followed by a ‘dry’ and a 
‘very dry’ year in 201 1 and 2012 (P = 379.2 mm in 201 1 and P = 288.6 mm in 2012). 

Details concerning the soil and vegetation conditions at the M-III site are provided in 
Wigneron et al. (2012). The field-site observed with ELBARA-II is typical of vineyards in the 
VAS region (the spacing between each plant is ~ 2 m and that between each row is ~ 3 m). 
Two field experiments in 2007 and 2010 led to similar values of the maximum Leaf Area 
Index, LAImax ~ 2.2. To monitor the time variations in the vegetation characteristics over the 
growing season, we used the NDVI index from the MODIS products (16 day NDVI 
composite of 250 m MODIS data; MODIS (2010)). As the field was large enough (larger than 
300 m x 300 m), it can considered that the MODIS NDVI time variations are representative of 
the vegetation conditions as seen by the ELBARA-II radiometer operated at the M-III site. 

In order to monitor the evolution of the surface roughness over time, field measurements were 
made by means of measuring mechanically two-dimensional profiles of the ground surface. 
For this purpose, a 2 m needle board with 201 needles, movable in the vertical direction and 
with 1 cm spacing between needles was used (Mialon et al. (2012)). The needle board was 
leveled and placed on the ground such that the needles were allowed to fall until they touched 
the soil surface. Subsequently, photos of the profile created by the needle heights were taken 
and digitized to compute soil roughness parameters. On each of the seven days during 2012 
when roughness measurements were performed, approximately 8 to 12 profdes were taken 
within the ELBARA-II footprints. Different locations and orientations (perpendicular and 
parallel to the vegetation rows) were considered in computing representative information on 
the standard deviation of soil surface height (So, cm), and correlation length (Lc, cm). Time 
variations in the average values of Sd and L c are shown in Fig. 1 for the seven days of 
measurements in 2012. The corresponding annual mean values are (Sd) = 2.2 cm, (Lc) 
= 6.2 cm. 


2.3 L-MEB modelling and inversion 
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The data set considered as a reference in this study was obtained using the 2-P L-MEB 
inversion approach to obtain retrievals of SM and xnad (Wigneron et al., 2000). There are 
many reasons to use this retrieved data set as a reference. 

First, the SM data set retrieved from tower-based remote sensing observations can be 
considered as representative of the SM conditions over the whole ELBARA-II footprint (this 
is usually a complex task using field probes distributed within the field). Second, the 2-P L- 
MEB method, based on multi-angular observations, has been validated in many studies 
against experimental data sets for a variety of soil and vegetation conditions (Wigneron et ah, 
1995, 2007; Parde et ah, 2003, 2004, Saleh et ah, 2006; Panciera et ah, 2009; Cano et ah, 
2010; Schlenz et ah, 2012, etc.), and its accuracy and robustness has been evaluated 
theoretically (Wigneron et ah, 2000). The 2-P L-MEB method is currently implemented in the 
official SMOS retrieval algorithm (Kerr et al, 2012). Third, the 2-P L-MEB approach has the 
advantage of providing retrievals of optical depth at nadir (tnad). These latter values were 
used to calibrate the relationships between tnad and NDVI, which are required in the SCA-E1 
and SCA-V algorithms. Moreover, it can not be considered that one method can benefit from 
the use of 2-P L-MEB retrieval method as a reference: the equations of the L-MEB model, 
used in the 2-P L-MEB approach, are also the basis of the SCA-E1, SCA-V, DCA and LPRM 
algorithms. 

A detailed description of the L-MEB model is given in Wigneron et al. (2007) and a brief 
summary of the main L-MEB equations and of additional parameterizations developed since 
2007 is given in the following. The L-MEB model is based on a zero-order solution of the 
radiative transfer equations: the so called x-co model, where the optical depth x accounts for 
extinction effects within the canopy and the effective scattering albedo co (-) accounts for 
scattering effects (Mo et al., 1982; Kurum et al., 2013). To incorporate the SMOS multi- 
angular feasibility, several additional parameterizations are used in L-MEB to account for 
effects of the vegetation structure and soil roughness on L-band brightness temperatures 
emitted from vegetated land surfaces. 

In local thermal equilibrium the emissivity eop of the ground at horizontal (p = H) and vertical 
(p = V) polarization is related to the corresponding reflectivity rop of the soil (the ground) 
observed at the angle 0: 


eop (6) - 1 - Asp (0) 
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The soil reflectivity r GP can be expressed as the reflectivity r*op of a specular surface and the 
roughness model parameters Qr, Hr and Nrp as: 

r G p (6) = [(1 -Qr) r* GP (0) + Q R r* GQ (0)] exp (- H R cos NRP (0)) (2) 

In this equation, Hr parameterizes the intensity of the roughness effects, Qr parameterizes the 
polarization mixing effects, and Nrp is used to account for the specific effects of roughness on 
the trend of soil reflectivity r G p as a function of incidence angle and polarization. The 
reflectivity of a specular surface r* GP was computed using the Fresnel equations as a function 
of 0 and of the effective soil dielectric permittivity e G . The latter was computed from soil 
moisture SM, soil effective temperature T G , and from the clay fraction using the dielectric 
mixing model of Mironov et al. (2012), referred to as the ‘Mironov’ model in the following. 
This is in contrast to the earlier study Wigneron et al. (2007), where the Dobson model 
(Dobson et al., 1985) was used to estimate e G . 

We used the recent results of Lawrence et al. (2013) to estimate the values of the roughness 
model parameters (Qr, Hr and Nrp). These parameters were assumed as constants in time, 
and therefore computed from the annual average value (Sn) of the standard deviation of the 
soil surface height and the corresponding annual mean (L c ) of the correlation length (Fig. 1). 
To be consistent with the general approach considered for SMAP we assumed that Nrv = Nrh 
= 0 (O’Neill et al., 2013). On that assumption, the roughness parameters Hr and Qr were 
computed as (Lawrence et al, 2013): 

Hr = 1.762 (1 - exp (-Z s /1.85)) and Q R = 0.05 Hr (3) 

where Z s =(Sd) 2 / L c (cm) 

Considering the annual mean values (So) = 2.2 cm and (Lc) = 6.2 cm measured over the M-III 
site in 2012, we obtained Z s = 0.78 cm, H R = 0.606, Q R = 0.0303. 

In this study, we considered a composite soil-vegetation surface temperature Tqq for the 
effective temperature Tq of the ground (the soil) and the vegetation canopy 7) . The composite 
effective temperature Toe of the ELBARA-II footprints was computed from the ERA- 
INTERIM 0-7 cm soil temperature product (Te- 07 ). ERA-INTERIM is the latest ECMWF 
(European Centre for Medium-Range Weather Forecasts) global atmospheric reanalysis of the 
period 1979 to the present (Dee et al., 2011) with a temporal resolution of 3 hours and a 
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282 spatial resolution of 0.75° (corresponding to about 100 km resolution over the VAS site). The 

283 accuracy of this estimate was considered to be sufficient in several studies investigating SM 

284 retrievals from L-band observations (Parde et al., 2004; Wigneron et al., 2012). 

285 As noted above, we used the x - co model to compute the upwelling emission (T B ) from the 

286 two layer soil-vegetation medium. T B p (p = H, V) is the sum of three terms: (1) the direct 

287 upwelling vegetation emission, (2) the downwelling vegetation emission reflected by the soil 

288 and attenuated by the canopy layer, and (3) upwelling soil emission attenuated by the canopy: 

289 T B p = (l-<Dp) (1-yp) (1 + y P rep) Tc + (1-rop) yp To (4) 

290 where Tg = T c = Tgc = T e _ 0 7 is assumed in this study, and r G p is the soil reflectivity computed 

291 with (2) and (3). y P is the vegetation attenuation factor which is related to the optical depth ip 

292 as (Beer’s law): 

293 y P = exp( - xp / cos 0) (5) 

294 To account for vegetation anisotropies, the optical depth xp(0) at the observation angle 6 is 

295 expressed with a parameterization involving the optical depth xnad at nadir (0 = 0°) : 

296 xp(0) = Xnad (sin 2 (0).tt P + cos 2 (0)) (at p = V, H) (6) 

297 The parameters tty (-) and tt B (-) account for the angular dependence of xp(0). As found in 

298 Wigneron et al. (2012), we considered that tt B = 1 (default L-MEB value) and that the tty 

299 parameter is free in the retrieval process, to account for the effects of the vine stocks, with a 

300 preferential vertical orientation. So in reality, a 3 -Parameter retrieval approach is made in this 

301 study, but the notation 2-P is kept, as only SM and x can be considered as variables of interest 

302 for applications. 

303 The values of the effective scattering albedo cop were found to be close to zero over most of 

304 the non-forested vegetation covers (Grant et al., 2008; Kurum et al., 2013). The value of cop 

305 was set equal here to 0.02 for both polarizations. A summary of the values of the soil and 

306 vegetation L-MEB parameters used in this study over the M-III site and described above is 

307 given in Table 1. 

308 The 2-P L-MEB inversions were based on bi-polarization and multiangular T B measurements 

309 using a minimization procedure of a cost function evaluating the difference between the L- 
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MEB simulations and the Tb measurements (Wigneron et al., 2000, 2007, 2012). The 
retrievals were based on ELBARA-II Tb data acquired with the automated elevation scans 
(section 2.1) performed for the observation angles 6 = 30°, 35°, 40°, 45°, 50° (corresponding 
roughly to the limit of validity of L-MEB at large incidence angles). As noted above, only Tb 
measurements made at 6 am and 6 pm will be considered in this study. Especially for the 
measurements at 6 am temperature gradients across the vegetation and the soil are minimal 
(Kerr et al., 2001). 


2.4 Description of the different SM retrieval methods 

As mentioned in the introduction, seven SM retrieval approaches were evaluated and 
compared in this study: the four methods considered presently in the SMAP ATBD for the 
passive-only product and three regression methods (described in Saleh et al (2006) and 
Mattar et al. (2012)) developed in the context of SMOS. The retrieved SM values were 
compared to a 'reference' SM data set obtained from the 2-P L-MEB inversion, which was 
assumed to be representative of the SM values over the ELBARA-II footprint. The seven SM 
retrieval approaches are described in the following sections. As is the case for the 2-P L-MEB 
method, these seven methods use the x-co radiative transfer model (described above) to 
account for the vegetation effects and they all assume tnad is independent of polarization and 
incidence angle (x v (0°) = Xh(0° ) = x N ad)- They are based on the same equation (1) to model 
the roughness effects, considering that Nrv = Nrh = 0. Furthermore, as implemented here they 
all use the ‘Mironov’ equations to compute the effective soil dielectric permittivity £g. All of 
the parameters listed in Table 1 for the 2-P L-MEB method are accounted for in the seven SM 
retrieval methods considered. Only a very brief description of the SCA-H, SCA-V, DCA and 
LPRM methods will be given here as a detailed description of these methods is available in 
the initial release of the ATBD. All these four methods were applied to the ELBARA-II T B 
data at the incidence angle of 40° corresponding to the SMAP observations. A summary of 
the input variables required for the seven different retrieval methods, as well as for the 
reference algorithm 2-P L-MEB, is given in Table 2. 


2.4.1 Single Channel Algorithms (SCA-H and SCA-V). 
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The Single Channel Algorithm (SCA-H), based on horizontally polarized Tb observations, is 
the current SMAP baseline, but the same algorithm can also be applied to vertically polarized 
Tb data (SCA-V). In SCA-H, brightness temperatures are converted to emissivity using a 
surrogate for the temperature of the emitting surface layer (in this study, the soil temperature 
provided by ECMWF (T E _ 0 7 ) is used). The derived emissivity is corrected for vegetation and 
surface roughness to obtain the soil emissivity. Finally, a dielectric mixing model (the 
‘Mironov’ model in this study) is used to obtain soil moisture SM from the soil dielectric 
constant Sg using the Fresnel equations. 

In this investigation, SCA-H and SCA-V are based on the same corrections of vegetation 
(using the x-co model), and soil roughness effects (using the H R and Q R parameters) as those 
used for the 2-P F-MEB method. 

t\ad is estimated from the vegetation water content (VWC) as 

x NAD = b . VWC (7) 

where b is a proportionality factor mainly depending on the vegetation structure. 

For SMAP, values of b will be provided by means of a land cover look up table and the 
baseline approach utilizes a set of land cover-based equations to estimate VWC from values 
of NDVI. The following equation is used for cropland (O’Neill et al, 2013): 

VWC = (1.9134 x NDVI 2 - 0.3215 x NDVI) + Stemfactor x (NDVI ref - 0.1) / (1 - 0.1) 

( 8 ) 

where Stemfactor parameter is the product of the average height of a land cover class and the 
ratio of sapwood area to leaf area; NDVI re f is assumed to be equal to the maximum value of 
NDVI time series (the value of NDVI re f was set equal here to 0.4696 from the analysis of the 
MODIS NDVI observations over the 2010-2012 period). In this study, the b and Stemfactor 
parameters were calibrated prior to the inversion process, as described in Section 2.5. 

2.4.2 The Dual Channel Algorithm (DCA) 

The Dual Channel Algorithm (DCA) is an extension of the SCA and uses both H-polarized 

and V-polarized Tb observations to simultaneously retrieve SM and VWC (O’Neill et al, 

2013). As in the 2-P F-MEB algorithm, the SM and tnad variables are adjusted iteratively 

until the root mean square difference between the simulated and observed Tb is minimized. 

There are differences between 2-P F-MEB and DCA algorithms. Firstly, T B data at 0 = 40° 

are used for DCA, while multiangular data are used for 2-P F-MEB. Secondly, the tty 
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parameter (accounting for an angular dependence of x) is retrieved in 2-P L-MEB, while DCA 
does not account for this dependence. Except for the tt v and tt H parameters, all vegetation and 
soil parameters used in DCA are the same as those used in the 2-P L-MEB method (Table 1). 

2.4.3 Land Parameter Retrieval Model (LPRM) 

The LPRM approach uses an analytical solution for the derivation of the vegetation optical 
depth. This solution uses the Microwave Polarization Difference Index (MPDI) and the 
observed surface emissivity (en and ey) as input and is based on the assumption that the 
values of the vegetation optical depth are the same for both polarization (xy = xh). The MPDI 
index is calculated from the brightness temperature at El- and V polarizations as follows 
(Meesters et al., 2005): 


MPDI = (T bv - T bh ) / (T B y + T bh ) (8) 

Then based on equation (4) of the x-co omega model, soil moisture is retrieved using an 
optimization routine that minimizes the RMSE between the modelled and observed El- 
polarized brightness temperatures. As for SMOS, the vegetation optical depth at this 
optimized soil moisture value is an additional retrieval result. As noted in O’Neill et al. 
(2013), the LPRM was implemented on multifrequency satellites such as AMSR-E, where 
also the Ka-band V-polarized chamiel is used to retrieve physical temperatures of the scene 
observed. This latter can also be estimated from re-analysis or near real time data from 
weather prediction centres (Parinussa et al., 2011), as is done in the current SMOS SM 
retrieval algorithm (Kerr et al., 2012). Only a few studies (e.g. de Jeu et al., 2009) have 
examined the applicability of this model at L-band frequencies, although the analysis of 
SMOS data with LPRM is currently underway. All detailed equations of the LPRM approach 
are given in (Owe et al., 2001; Meesters et al., 2005, Owe et al., 2008, de Jeu et al., 2009, 
Chung et al., 2013). As for DCA, except for the tty and ttn parameters which are not relevant 
here, all vegetation and soil parameters used in LPRM are the same as those used in the 2-P 
L-MEB method (Table 1). 

2.4.4 Linear regression methods (Saleh et al., 2006; Mattar et al., 2012) 
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Two methods based on regression equations developed by Saleh et al. (2006) and Mattar et al. 
(2012) were evaluated in this study. Both methods were numerically derived from the 
equations of the x-co model assuming, as for LPRM, that the value of the effective scattering 
albedo is cop = 0, and that the values of optical depth xp are the same for both polarizations p = 
H, V. These methods are physically-based. However, as the development of an analytical 
fonnulation would be complex, most of the time they are used as regressions methods. As 
shown by Saleh et al. (2006), a key interest in these regression methods is that they can be 
used for varying roughness and vegetation conditions over time: no additional information 
about temporal changes in these two state variables (such as NDVI or LAI for vegetation for 
instance) is required. These regression methods have been used in several studies based on in 
situ, airborne or spaceborne (SMOS) observations (Albergel et al., 201 1; Parrens et al., 2012; 
Calvet et al., 2011, etc.) 

The method of Saleh et al. (2006) can be applied to observations made either at the two 
incidence angles 0i and 02 (referred to as ‘Saleh’ bi-angular): 
ln(SM) = a 2 ln(T P (0i)) + ai ln(T P (0 2 )) + a 0 (0i, 0 2 , p) 

(9) 

or to bi-polarization observations made at one observation angle 0 (referred to as ‘Saleh’ bi- 
polarization): 

ln(SM) = b 2 ln(T H (0)) + b, ln(T v (0)) + b 0 (0) 

( 10 ) 

where Tp(0) is the reflectivity of the soil-vegetation system at polarization p (p=V or p=H), 
defined as 

T P (0)= 1 -Tbp(0)/T gc (11) 

where the composite soil vegetation surface temperature Toe was estimated from the ERA- 
INTERIM 0-7cm soil temperature product (Te- 07 ). 

The method of Mattar et al. (2012) is very similar and can be written as (referred to as 
‘Mattar’): 

In (SM) = c 2 In (r P (0)) + ci NDVI + c 0 (0, p) 

( 12 ) 

where the NDVI is considered here as a proxy for optical depth, as in the SCA-H and SCA-V 
methods. 
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In the above equations (9), (10) and (12), the parameters (a 0 , ai, a?), (bo, bi, b2) and (co, Ci, C2) 
are regression coefficients, which are assumed to be constant in time and have to be calibrated 
over each pixel. In this study, in the ‘Saleh bi-polarization’ equation (10), we used the 
observation angle 6 = 40° as used in the other retrieval methods. In the ‘Saleh bi-angular’ 
equation (9), we used H-polarized bi-angular observations at 0| = 30° and 0 2 =50°. In the 
‘Mattar’ equation (12), we used H-polarized observations at 0 = 40°. These latter 
configurations were found to be the best for SM retrievals (results not shown here). 


2.5 Method calibration 

In this study, the SCA-V, SCA-H, DCA and LPRM methods were based on the L-MEB 
model parameters given in Table 1. In addition, some model parameters specific to some 
methods had to be calibrated. The DCA and LPRM methods did not require any additional 
calibration. Conversely, in the SCA-V and SCA-H methods, the two parameters b and 
Stemfactor, used to link NDVI and optical depth, had to be calibrated. Moreover, the three 
‘regression’ methods ‘Saleh bi-angular’, ‘Saleh bi-polarization’ and ‘Mattar’ did not require 
any L-MEB parameters but required the calibration of three coefficients (a,), (bl) or (c,) (i = 0, 
1 and 2) used in equations (8), (9) and (11), respectively. 

The calibration of the above parameters and coefficients was perfonned three times, using one 
year of data for calibration and the two other years for validation. To calibrate the b and 
Stemfactor parameters in SCA-H and SCA-V, a multilinear regression method was used to fit 
the optical depth derived from equations (6) and (7) to the ‘reference’ optical depth t N ad 
retrieved from the 2-P L-MEB method. The obtained values for all three calibration years 
(2010, 201 1 and 2012) are given in Table 3. 

Similarly, to calibrate the three coefficients in the regression equations of the ‘Saleh bi- 
angular’, ‘Saleh bi-polarization’ and ‘Mattar’ methods, a multilinear regression method was 
used to minimize the difference between the retrieved SM derived from equations (9), (10) or 
(12) to the ‘reference’ SM values retrieved from the 2-P L-MEB method. The obtained values 
of the coefficients for all three methods and all three calibration years (2010, 201 1 and 2012) 
are given in Table 3. 


3. Results 
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3.1. Reference values of SM and tnad 


As outlined above, the ‘reference’ values of soil moisture (SM) and optical depth at nadir 
(tnad) were retrieved from the multiangular Tb data measured by the ELBARA-II instrument. 
The Tb measured 0 — 40° for the time period 2010-2012 are shown in Fig. 2. A clear seasonal 
cycle in the Tb time-series can be seen, with maximum values of Tb during summer and lower 
Tb values during winter. This annual cycle is related to the vegetation growth cycle, 
beginning in April and ending in November, and to the soil moisture conditions, which are 
generally drier during the summer period. 

However, as already noted in Section 2.2, significantly wetter/drier conditions were 
encountered in 2010/2012, respectively, which is reflected in the observed Tb trends over the 
MELBEX-II site with lower values during summer 2010 compared to summer 2012. Based 
on these T B observations, the retrieved values of SM and t NA d were computed from the 2-P L- 
MEB method and they are illustrated in Fig. 3a-b. As discussed in Jackson et al. (2012), 
conditions of standing water during or shortly after intensive rainfalls should be flagged. In 
this study, to avoid these conditions, all retrieved values of SM which were found to be larger 
than the saturation value SMsat were not considered (SMsat was set equal to 0.5 m 3 /m 3 over 
the M-III site as computed by Juglea et al. (2010)). Note that due to this data fdtering, the 
number of SM data used in the comparison may vary slightly from one approach to the other. 

In accordance with the above-discussed Tb trends one can see that rainy conditions led 
generally to higher values of SM throughout the year in 2010 and during the winter period in 
2011 and 2012 (Fig. 3a). Drier conditions during the second half of 2011 and 2012 led to 
rather long time intervals of lower SM values. 

The vegetation cycle could be clearly distinguished from the time variations in both the 
optical depth at nadir (tnad) and NDVI index obtained over the 250 m MODIS pixel 
including the M-III vineyard (Fig. 3b). Relatively similar maximum values of tnad were 
retrieved during the summer of all three years (maximum values of tnad are close to 0.24 in 
2010 and close to 0.22 in 201 1 and 2012). During the winter period, after vine pruning and 
defoliation, values of tnad close to 0.05 were retrieved for all three years. This latter value 
corresponds to the estimated value of the optical depth (x stock) of vine stocks (Schwank et 
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al., 2012; Wigneron et al., 2012). Superimposed on the long term trend of xnad, short-time 
changes in the time variations of tnad can be noted. It is likely that these apparent fluctuations 
result from unaccounted for changes in the roughness conditions over the field as discussed in 
Patton and Hornbuckle (2013) and Jackson et al. (2012) for SMOS observations. It can be 
noted too that very low values of xnad were retrieved during a short period of time in May of 

201 1 and 2012, just before the vine vegetation growth. We assumed that this could be caused 
by specific effects during this period related to soil roughness or to vegetation structure. For 
instance, this effect could be linked to lower roughness conditions in relation to field works in 
May. As the roughness parameterization is set as constant over the 3 year period, actual lower 
roughness conditions in the field would lead to retrievals of lower x N ad values and, to a lesser 
extent, higher SM values. Our field observations of roughness for the year 2012 (Fig. 1) are 
not accurate enough to confirm clearly this assumption but they seem to be leaning in that 
direction. 

A maximum value of NDVI is reached in the middle of July (~ DoY 200): NDVImax ~ 0.45 
in 2010 and 2011 and NDVImax ~ 0.36 in 2012. It is likely the lower value of NDVImax in 

2012 can be related to the drier conditions during that year. In comparison with the year 2011, 
it seems that the very dry conditions during 2012 impact the NDVI values, but do not impact 
the time variations of Xnad considerably. 

A scatter plot of the retrieved values of the optical depth xnad versus the NDVI index is 
shown in Fig. 4. It can be seen that the results are generally consistent from one year to the 
other. One specific pattern can be noted in 2011; it corresponds to very low values of xnad 
retrieved while vegetation is fully developed (NDVI « 0.45), which was already discussed 
above. 


3.3 Comparison of SM Retrievals 

The retrieved values of SM from all retrieval methods presented in section 2 were compared 
to the reference SM values retrieved with the 2-P L-MEB method applied to the 
measurements performed during the years 2010-2012. A summary of this comparison is given 
in Table 4, in terms of coefficient of determination (R ), bias (m / m ), RMSE (m / m ) and 
unbiased RMSE (ubRMSE, m 3 /m 3 ) as defined by Entekhabi et al. (2010). To illustrate the 
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results, scatter plots of retrieved SM values versus ‘reference’ SM values are given for all 
methods considered in this study (Fig. 5). 

All five methods requiring a calibration step, e.g. SCA-V, SCA-H, ‘Saleh’ bi-angular, ‘Saleh’ 
bi-polarization and ‘Mattar’ (the calibration was made using one year and the evaluation with 
the two other years), provided SM retrievals that were in good agreement with the ‘reference’ 
SM data (R 2 is generally higher than 0.90, and the RMSE is lower than 0.045 m 3 /m 3 ). If we 
consider the years used for calibration, best performances in terms of R 2 for all four methods 
were obtained when year 2010 (corresponding to rather ‘wet’ conditions) was used for 
calibration, while lower performances were obtained using the year 2012 (corresponding to 
‘very dry’ conditions) for calibration. Results for the year 2011 are generally close to those 
obtained for the year 2010. A closer inspection shows that both the SCA-V and the SCA-H 
methods provide generally very similar performances in SM retrievals (the SCA-V method 
providing a slightly better accuracy in terms of R 2 , bias, RMSE and ubRMSE). The three 
methods based on regression equations (‘Saleh’ bi-polarization, ‘Saleh’ bi-angular, and 
‘Mattar’) provided very similar results too. Slightly lower performances were obtained for the 
‘Mattar’ method (especially when using the year 2010 for calibration), while best 
performances were obtained for ‘Saleh bi-angular’. Considering the ubRMSE criteria, the 
performances of the SCA and the regression methods were even closer. Except for the DCA 
and LPRM algorithm, the ubRMSE is always around or below the target accuracy for SMAP 
of 0.04 m 3 /m 3 . This is a direct consequence of the fact that values of the bias were found to be 
higher for the SCA methods (bias ~ 0.020 m 3 /m 3 ) than for the regression methods (bias ~ 
0.010 m 3 /m 3 ). 

As could be expected, results obtained from methods which did not require parameter 
calibration (DCA and LPRM) provided results with a lower accuracy: the RMSE was similar 
for both methods (RMSE ~ 0.55 m 3 /m 3 ), while slightly better R 2 values were obtained for 
DCA (R 2 = 0.79) than for LPRM (R 2 = 0.725). For both methods, the bias in the retrievals 
was found to be very low (bias = 0.021 m 3 /m 3 for DCA, and bias = 0.013 m 3 /m 3 for LPRM). 

The scatter plots (Fig. 5) showing the comparison between retrieved SM values versus 
‘reference’ SM values are given to illustrate these different results. For methods requiring 
calibration (SCA-V, SCA-H, ‘Saleh’ bi-angular, ‘Saleh’ bi-polarization and ‘Mattar’), we 
used the year 2010 in Fig. 5 (this year provided best performances in terms of R 2 ). Note that 
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the number of data used in the comparison may vary from one approach to the other. This can 
be explained by two reasons. First, for DCA and LPRM, the comparison was made over three 
years (2010 - 2012), while it was made over two years (2011 -2012) for the other methods. 
Second, retrieved SM values larger than the saturation value SM SA t (SM S at = 0.5 m 3 /m 3 ) 
were removed in the comparison (a very low number of observations was concerned by this 
filtering). 

It can be seen that a very low bias was obtained generally. Flowever, in wet conditions, the 
methods LPRM and DCA provided underestimated SM values (for SM > 0.3 m 3 /m 3 ); while 
the ‘Saleh’ and ‘Mattar’ methods provided overestimated SM values (for SM > 0.2 m 3 /m 3 ) 
with respect to the reference SM. For DCA and LPRM (methods with do not require any 
calibration), it can be seen that the SM retrieval performances are lower in a small SM 
interval, for values of SM comprised between ~ 0.1 and 0.15 m 3 /m 3 . These SM conditions 
generally correspond to periods of vegetation growth at the end of spring and of full 
vegetation development in the summer period. 


4. Discussion and conclusion 

This study presents an inter-comparison of several SM retrieval methods based on a three year 
data set of passive L-band microwave observations acquired over a vineyard site at the VAS 
site. 

A careful interpretation of the results should be made, and the results cannot be easily 

generalized to operational applications for spaceborne sensors. We will discuss these different 

aspects and the main conclusions of the study in the following. First, it is important to 

consider that the results were obtained at the field scale and over only one type of vegetation 

(a vineyard canopy) with some specific features (no litter layer, relatively low LAI and 

biomass conditions, frequent agricultural practices leading to changes in soil roughness, etc.). 

Several effects related to changes in the soil roughness conditions or in the vegetation 

structure (in relation with the crop growth and the agricultural practices) may have a 

significant impact on the results of the present study. It is likely that the impact of these 

effects would average out and, therefore, become much less important if we had considered 

larger footprints of spaceborne radiometric observations, including a large variety in the types 

of vegetation (natural or cultivated canopies), in the soil conditions and in the agricultural 
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practices. For instance, specific effects related to the vegetation structure could be revealed 
over the vineyard field and the values of optical depth for both polarizations (xh(0) and xv(0)) 
could not be considered as equal for that canopy type (Wigneron et al., 2012). This result has 
frequently been obtained from in situ radiometric observations (Parde et al., 2003, 2004; 
Wigneron et al., 2004) but it has never been noted, to our knowledge, from spaceborne 
observations. For instance, Owe et al. (2001) found that xy = xh over test sites in the US over 
a variety of land covers based on SMMR (Scanning Multichannel Microwave Radiometer) 
observations at C-band. It is likely that these vegetation structure effects can be a limitation 
for presented evaluation of the methods, which all assume xh = xy. So, several similar studies 
based on in situ observations over a variety of vegetation types are required to provide a more 
in-depth evaluation of the method performances. 

It should be noted too that the performances of the different methods cannot be compared 
directly as some methods had to be calibrated while some methods did not require any 
parameter calibration step (DCA and LRPM). The two methods SCA-H and SCA-V, require 
the calibration of the relationship between optical depth and a remotely sensed vegetation 
index (NDVI); the three methods based on regression equations, ‘Saleh bi-angular’, ‘Saleh bi- 
polarization’ and ‘Mattar’, require the calibration of three coefficients. This calibration step 
could be done in the present study as we considered that a ‘reference’ data set describing the 
time variations in SM and xnad (and derived from multi-angular observations) was available 
from the ELBARA-II tower-based observations. However, for operational spatial 
applications, it is generally very difficult to obtain such a reference data set. 

In spite of the limitations discussed above, some key results obtained in this study from 
tower-based observations could be of value to future operational applications. It was found 
that the two methods, which did not require any a priori calibration (DCA and LPRM) could 
provide good SM retrievals and have relatively similar performances (R 2 ~ 0.72-0.79; 
RMSE ~ 0.054-0.58 m 3 /m 3 ) over the three year period. The methods requiring parameter 
calibration (two parameters in SCA-H and SCA-V; three coefficients in the three regression 
methods) provided results closer to the reference: for instance the R 2 coefficient increased 
generally to values larger than 0.90 for all methods. The methods which require additional 
information concerning the vegetation development (the NDVI variable is required in the 
SCA-H, SCA-V and ‘Mattar’ algorithms) provided slightly lower performances when year 
2012 was used for calibration. For that year the NDVI values were lower than for the two 
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other years (maximum NDVI values ~ 0.45 in 2010 and 201 1 and ~ 0.36 in 2012), while the 
maximum values of tnad were found to be relatively similar over all three years ( ~ 0.22 - 
0.24). It is likely that nonlinearities between xnad and NDVI led to these slightly lower 
performances in SM retrievals for the year 2012 for the SCA and ‘Mattar’ algorithms. 

In the present study, the computed performances are “optimal” performances as it is assumed 
that a good parameter calibration can be made from a SM data set which can be considered as 
a reference. This calibration step was possible in this study based on in situ tower-based 
observations obtained over a homogeneous vineyard field, but this step is much more 
complex for operational applications based on space borne sensors. Several options are 
possible to calibrate these different retrieval methods for spaceborne applications. For 
instance, the reference SM or tnad values which are required in the calibration step can be 
estimated: 

(i) from networks of in situ measurement sites such as SCAN in the USA (Schaefer et al., 
2007), OZNET in Australia (Smith et al., 2012) or SMOSMANIA in France (Albergel et al., 
2012), etc. Then, based on results obtained over a variety of soil and vegetation conditions, a 
look up table providing the calibrated parameters as function of the land cover types can be 
built. 

(ii) from model re-analyses (ERA-Interim (Dee et al., 201 1) or MERRA Land (Reichle et al., 
2012) for instance), in regions where the simulated SM values can be considered to be 
accurate. As mentioned above, in a second step, a look up table can be built for a variety of 
land covers. 

(iii) by combining observations from different remote sensing sensors. For instance, the 
estimation of optical depth tnad retrieved from SMOS or other satellites (e.g. AMSR-2) could 
be used to calibrate the vegetation parameters required in the SCA-H and SCA-V algorithms 
(Lawrence et al., 2014). 

Future work will consider these different options to evaluate the retrieval capabilities of the 
different methods requiring calibration (SCA, ‘Saleh’ or Mattar’) for operational applications 
based on spaceborne sensors. 
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842 Table 1. 

843 L-MEB soil and vegetation parameters over the M-III vineyard (VAS site). All these parameters, 

844 except tt H and tt v which are specific to L-MEB, are valid for the other SM retrieval methods. 



Unit 

Value or used Model 

Soil dielectric 

(-) 

Mironov et al. (2012) 

permittivity (e G ) 
Clay fraction 

(-) 

0.26 (in situ measurements; Juglea et al., 2010) 

H 

a 

II 

H 

n 

II 

H 

o 

o 

K 

ECMWF ERA Interim temperature (T E _ 0 7 ) 

Hr 

(-) 

0.6060 (calibrated, Lawrence et al., 2013) 

Qr 

(-) 

0.0303 (calibrated, Lawrence et al., 2013) 

Nrh 

(-) 

0 

N rv 

(-) 

0 

tt H 

(-) 

1 

tty 

(-) 

Free parameter in the retrieval process 

CO 

(-) 

0.02 

X 

(-) 

Free parameter in the retrieval process 

SM 

m 3 / m 3 

Free parameter in the retrieval process 
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847 Table 2. 


848 Input variables required in the different retrieval algorithms 

849 


Algorithm 

Input variables 

SCA-H 

Tbh(9=40°) 

ECMWF temperature (T E _o 7 ) 
NDVI 

SCA-V 

Tbv(9=40°) 

ECMWF temperature (T E _o 7 ) 
NDVI 

DCA 

T bh (9=40°), T bv ( 9=40°) 
ECMWF temperature (T E _ 07 ) 

LPRM 

T bh (9=40°), T bv ( 9=40°) 
ECMWF temperature (T E _ 07 ) 

‘Saleh’ 

bi-polarization 

T B h(9=40°), T bv ( 9=40 o ) 
ECMWF temperature (T E _ 07 ) 

‘Saleh’, 

bi-angular 

Tbh(9=30°), T B h( 9=50°) 
ECMWF temperature (T E _ 07 ) 

‘Mattar’ 

T bh (9=40°) 

ECMWF temperature (T E _ 07 ) 
NDVI 
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852 


853 Table 3. 

854 Calibrated parameters of the different retrieval algorithms: one year (2010, 2011 or 2012) is used for 

855 calibration; the two other years are used for validation 


856 


857 


858 

859 


860 


SCA H/V, T B h( 0=4O) orT BV (0=4O) 


Calibration 

b 

Stemfactor 

2010 

0.61679 

0.20874 

2011 

0.31756 

0.44014 

2012 

0.92819 

0.05840 

‘Saleh bi-angular’, T BH (0=3O), T BH (0=5O) 


Calibration 

a 0 

ai 

a 2 

2010 

1.4171 

-0.3560 

0.8374 

2011 

1.0972 

-0.2806 

0.2613 

2012 

2.2857 

-1.5674 

0.1300 

‘Saleh bi-polarization’, T BH (0=4O), T BV (0=4O) 

Calibration 

bo 

bi 

b 2 

2010 

0.3524 

0.7734 

1.1401 

2011 

0.2595 

0.6208 

0.4879 

2012 

-0.3914 

1.1927 

0.7263 

‘Mattar’, T BH (0=4O) 

Calibration 

Co 

Cl 

c 2 

2010 

1.2530 

0.9491 

0.9147 

2011 

0.9844 

0.5748 

0.3702 

2012 

1.0954 

2.6578 

0.0183 


861 

862 


863 
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Table 4. 


865 Performances of the different SM retrieval algorithms in terms of coefficient of determination (R 2 ), 

866 bias (m 3 /m 3 ), RMSE (m 3 /m 3 ) and ubRMSE (m 3 /m 3 ). For SCA-H, SCA-V, ‘Saleh bi-angular’, ‘Saleh 

867 bi-polarization’ and ‘Mattar’, one year (2010, 201 1 or 2012) is used for calibration; the two others are 

868 used for validation. For LPRM and DCA, no calibration is required. 

869 


Method 

Calibration 

Validation 

R 2 

Bias (m 3 /m 3 ) RMSE 
(m 3 /m 3 ) 

ubRMSE 

(m 3 /m 3 ) 

SCA-H 

2010 

2011,2012 

0.915 

-0.025 

0.050 

0.043 


2011 

2010, 2012 

0.905 

-0.041 

0.054 

0.035 


2012 

2010, 2011 

0.852 

-0.020 

0.056 

0.052 

SCA-V 

2010 

2011,2012 

0.928 

-0.014 

0.035 

0.032 


2011 

2010, 2012 

0.919 

-0.024 

0.040 

0.032 


2012 

2010, 2011 

0.861 

-0.010 

0.045 

0.043 


DCA 0.789 0.021 0.054 0.050 


LPRM 0.725 0.013 0.058 0.056 


Saleh 

Bi-angular 

2010 

2011 

2012 

2011,2012 
2010, 2012 
2010, 2011 

0.950 

0.941 

0.934 

0.004 

0.007 

0.009 

0.037 

0.028 

0.036 

0.037 

0.027 

0.035 

Saleh 

2010 

2011,2012 

0.946 

0.010 

0.040 

0.039 

Bi- 

2011 

2010, 2012 

0.924 

-0.001 

0.031 

0.031 

polarization 

2012 

2010, 2011 

0.920 

0.004 

0.033 

0.033 


Mattar 

2010 

2011,2012 

0.946 

0.009 

0.041 

0.040 


2011 

2010, 2012 

0.927 

-0.001 

0.030 

0.030 


2012 

2010, 2011 

0.869 

0.017 

0.048 

0.045 
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877 

878 

879 

880 

881 

882 

883 

884 

885 

886 

887 
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889 
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Figure Captions 

Fig. 1 Temporal variations in the standard deviation of soil surface heights So and correlation 
length Lc estimated from measurements during seven days in 2012 performed at the M-III 
vineyard field. The annual mean values are (So) = 2.2 cm, (Lc) = 6.2 cm. 

Fig. 2. Time-series of measured ELBARA-II Tb over the M-III vineyard during three years 
(2010-2012) at H (‘o’) and V (‘x’) polarizations and at the observation angle 6 = 40°. The T B 
data are acquired ~ every 30 minutes but only data measured at 6 am are shown. Diurnal 
precipitation P is represented with vertical lines 

Fig. 3. Soil moisture SM (a) and optical depth tnad (b) retrieved with the multiangular 2-P L- 
MEB method applied to the measurements at the M-III site. The diurnal retrievals are shown 
for 6 am and 6 pm, respectively. These retrieved values are considered as a reference in this 
study. Diurnal precipitation is represented with vertical lines. In Fig 3b, the time-series of 
NDVI index obtained over the 250 m MODIS pixel including the M-III vineyard is shown. 

Fig. 4. Scatter plot of retrieved values of the optical depth tnad, retrieved with the 
multiangular 2-P L-MEB method, versus the NDVI index obtained over the 250m MODIS 
pixel including the M-III vineyard. Retrieved values of Tnad computed at 6 am and 6 pm are 
used. 

Fig. 5. Scatter plots of the retrieved SM values versus the reference SM values for all 
methods: SCA-H (a), SCA-V (b), DCA (c), LPRM (d), ‘Saleh’ bi-angular (e), ‘Saleh’ bi- 
polarization (f) and ‘Mattar’ (g). Retrieved values of SM are computed at 6 am and 6 pm. In 
Fig. 5a-b-e-f-g, retrieved values of SM for years 201 1 and 2012 are shown (the year 2010 was 
used for calibration). In Fig. 5c-d (for DCA and LPRM) retrieved values of SM for years 
2011, 2012 and 2013 are shown (no calibration was required). 
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895 Fig. 1 Temporal variations in the standard deviation of soil surface heights So and correlation 

896 length Lc estimated from measurements during seven days in 2012 performed at the M-III 

897 vineyard field. The annual mean values are (So) = 2.2 cm, (Lc) = 6.2 cm. 
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900 


Fig. 2. Time-series of measured ELBARA-II Tb over the M-III vineyard during three years 

901 (2010-2012) at H (‘o’) and V (‘x’) polarizations and at the observation angle 6 = 40°. The Tb 

902 data are acquired ~ every 30 minutes but only data measured at 6 am are shown. Diurnal 

903 precipitation P is represented with vertical lines. 
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905 Fig. 3. Soil moisture SM (a) and optical depth tnad (b) retrieved with the multiangular 2-P L- 

906 MEB method applied to the measurements at the M-III site. The diurnal retrievals are shown 

907 for 6 am and 6 pm, respectively. These retrieved values are considered as a reference in this 

908 study. Diurnal precipitation is represented with vertical lines. In Fig 3b, the time-series of 

909 NDVI index obtained over the 250 m MODIS pixel including the M-III vineyard is shown. 
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916 Fig. 4. Scatter plot of retrieved values of the optical depth tnad, retrieved with the 

917 multiangular 2-P L-MEB method, versus the NDVI index obtained over the 250m MODIS 

918 pixel including the M-III vineyard. Retrieved values of tnad computed at 6 am and 6 pm are 

919 used. 
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921 Fig. 5. Scatter plots of the retrieved SM values versus the reference SM values for all 

922 methods: SCA-H (a), SCA-V (b), DCA (c), LPRM (d), ‘Saleh’ bi-angular (e), ‘Saleh’ bi- 

923 polarization (f) and ‘Mattar’ (g). Retrieved values of SM are computed at 6 am and 6 pm. In 

924 Fig. 5a-b-e-f-g, retrieved values of SM for years 201 1 and 2012 are shown (the year 2010 was 

925 used for calibration). In Fig. 5c-d (for DCA and LPRM) retrieved values of SM for years 

926 201 1, 2012 and 2013 are shown (no calibration was required). 
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